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(57) ABSTRACT 

A color printing process, printing a color image in which 
out-of-gamut original colors are present. For each pixel 
defined by an original color which is determined to be out of 
gamut, a gamut remapping process is applied to map each 
pixel to a color which is within a printer gamut, remapping 
said pixels to colors within an output printer gamut. For a 
given set of gamut remapped pixels, gamut remapped pixel 
colors are compared with said original pixel colors, to derive 
a comparison metric. Using the comparison metric, a cor- 
rected set of gamut remapped colors is generated. The 
comparison metric may be subjected to an adaptive filtering 
process, which strengthens the comparison metric in high 
frequency image regions to increase its impact on the gamut 
remapped colors, and weakens the comparison metric in low 
frequency areas, to weaken its impact on the gamut 
remapped colors. 

16 Claims, 6 Drawing Sheets 




102 



FILTER 
SELECTOR 

'£ 

1 

ACTIVITY 
METRIC 



-302 



^400 



11/14/2003, EAST Version: 1.4.1 



U.S. Patent Nov. ll, 2003 Sheet 1 of 6 US 6,646,762 Bl 



i 


i Rl 1 IP 

TEXT 




LUMINANCE 






Y 




BLACK 

BACKGROUND 




SPATIAL POSITIONS 




FIG. 7 



LUMINANCE 
Y 1 



BLUE 
TEXT 



BLACK 

BACKGROUND 
TEXT 



SPATIAL POSITION X 

FIG. 2 



YELLOW 
STRIPE 



LUMINANCE 
Y 




SPATIAL POSITION X 

FIG. 3 



11/14/2003, EAST Version: 1.4.1 



U.S. Patent Nov. ll, 2003 Sheet 2 of 6 



US 6,646,762 Bl 



LUMINANCE 
V 



YELLOW 
STRIPE 

GREEN l l 
1 I RED 



SPATIAL POSITION X 

FIG. 4 



70 



20 



30 

I 



IAAAGE 




IPU 




PRINTER 


SOURCE 








GAMUT MAPPING 


"^24 





FIG. 5 



11/14/2003, EAST Version: 1.4.1 



U.S. Patent Nov. ll, 2003 Sheet 3 of 6 



US 6,646,762 Bl 




$2 



CM 

u 
u 



11/14/2003, EAST Version: 1.4.1 



U.S. Patent Nov. ll, 2003 Sheet 4 of 6 US 6,646,762 




FIG. 8 



11/14/2003, EAST Version: 1.4.1 



U.S. Patent Nov. ll, 2003 Sheet 5 of 6 




K 



FIG. 9 

11/14/2003, EAST Version: 



U.S. Patent Nov. ll, 2003 Sheet 6 of 6 



US 6,646,762 Bl 



COMPUTE 
SMALL AREA 

ACTIVITY a s 



300 




USE SMALL 
AREA FILTER 



YES 



COMPUTE 
LARGE AREA 

ACTIVITY a l 



306 



308 




YES 


USE SMALL 




AREA FILTER 



USE LARGE 
AREA FILTER 



372 



304 



310 



FIG. 10 



11/14/2003, EAST Version: 1.4.1 



US 6,646,762 Bl 

1 2 

GAMUT MAPPING PRESERVING LOCAL the luminance of the green side of the edge, while allowing 

LUMINANCE DIFFERENCES the yellow, which is within gamut, to print at its same 

luminance level, thus rendering it much more apparent and 

CROSS REFERENCE objectionable. 

Cross-reference is made to concurrently filed patent appli- 5 A similar problem occurs in antialiasing, where edge 

cation Ser. No. 09/435,256 entitled, "Gamut Mapping Pre- softening of saturated color edges results in printable colors, 

serving Local Luminance Differences with Adaptive Spatial and inconsistent out-of-gamut mapping of the saturated 

Filtering", by R. Balusubramanian, R. de Queiroz, R. Esch- colors causes the softened edge to stand out from the 

bach and Wengchen Wu. background and foreground colors. 

The present invention is directed to improving the selec- 10 U.S. Pat. No. 5,579,031 to Liang describes the use of a 
tion of printer colors where the original document selects printer model in the calibration, where the model is itera- 
colors that are out of gamut, and more particularly to an out tively updated to form a final LUT. While the use of a 
of gamut mapping method that preserves spatially local transform like the LUT described in itself inherently pro- 
luminance differences. vides color matching between the two devices considered, 

15 additional color gamut mapping may be used in the devel- 

BACKGROUND OF THE INVENTION opment of me LT j Tj by performing a preselected mathemati- 

Printing devices and display devices all inherently have a ca l manipulation that maps the input set of color values to he 

region of operation, sometimes referred to as a gamut. within the range of values from the printer model. Such 

Because such image reproduction devices are non-ideal, mathematical manipulation may be as simple as truncation 

they cannot reproduce every possible color that can be seen or as complicated as data compression, 

by a human. Gamut mapping is often used because the color U.S. Pat. No. 5,883,632 to Dillinger describes preserva- 

range that is possible to reproduce with a particular display tion of relative position by geometric transform on "hue 

is commonly not identical with the color range possible to pages". Device-space color specifications for two color 

reproduce with any selected printer. Thus, while both dis- 25 presentation devices are interrelated through the intermedi- 

play and printer have a large number of colors that are ary of at least one perceptual color space. This interrelation 

reproducible by both, there may be certain combinations of is performed by either tabulations or real-time software 

sets of color values that are outside the reproduction capa- processing, and in such a way that color, and color changes, 

bility of one or the other. However, images are composed for specified in relation to one device — in a control language of 

printing on display devices, and users expect printers to 3Q that device — are tracked in both the language and the actual 

reproduce images with their intent. Furthering the problem, performance of the other device. This mapping through 

different printing technologies and materials inherently pro- perceptual space can be used to produce the effect of 

vide different gamuts even among each class of devices. matching the two device gamuts to each other, so that the 

Gamut mapping serves to map pixels defined by colors full gamuts of both are in effect merged. When so used, this 

not printable by the printer or resulting from initial image 35 new kind of mapping preserves relative positions in color 

processing, into colors printable by a real printer. In doing s P ace — and tDUS tne capability of the tracking device to 

so, out-of-gamut, unprintable colors are mapped to printable maintain distinctions between colors seen on the source/ 

colors in accordance with some scheme designed to main- m P ut device that appear very closely similar. Such discrimi- 

tain color intent and aesthetic appearance. One of the most nation between adjacent colors can be supported essentially 

common schemes is to map any out-of-gamut pixel to its 40 throughout the gamut of both devices and in particular even 

closest in-gamut neighbor. Obviously this is not satisfactory, very near the lightness extrema and maximum-saturation 

because in a region of slight color variation, that is out- point °n a hue page. 

gamut, a number of similar colors may be mapped to the An article by T. Morimoto et al, "Objects* Color Match- 
same in-gamut color. Some textures, hi^ilights and shadows ings Based on Image Clustering" (Japan Hardcopy '98 p. 
in an original image will be lost. 45 371 ff), teaches image segmentation and color mapping for 

This problem leads to a generalization that a fundamental the different objects, 

problem of traditional gamut mapping processes is that they U.S. Pat. No. 5,903,275 to Guay describes a gamut 

are pixelwise operations that do not take neighborhood mapping in which saturated colors (i.e., out of gamut colors) 

effects into account. Looking at one example, illustrated at are mapped to device saturated values (device printable 

FIG. 1, when blue text is placed against a black background, 50 saturated colors). Colors inside the gamut are mapped to 

the original monitor or display has no difficult in providing "closest" and other colors to a blend between closest values 

an image that clearly distinguishes the blue and black calorimetrically and saturation. First, the printable colors are 

regions. However, the highly saturated blue color may be out measured is some device independent color space, such as 

of gamut of some printers, and black is commonly mapped L*a*b*. Next, the range of colors from some reference color 

to the printer's blackest reproduction. For large area 55 CRT, consisting of the colors from black to monitor 

reproduction, this maybe entirely acceptable. However, with saturates, and the colors between, then from saturates to 

the colors juxtaposed, looking at a luminance profile as in white, and the colors between, are calculated. The saturates 

FIG. 2, the luminance difference is decreased dramatically, are assigned colors based on what printable color a user 

rendering the appearance difference on the printed page would like to get when the user asks for a saturated CRT 

problematic. 60 color. This gives the printer's representation of the outer 

In another example, shown in FIG. 3, a saturated red/ range of the color gamut of the CRT. Next, a table is 

green edge that has been processed with JPEG compression generated that subsamples all colors available to the speci- 

and decompression, results in an intermediate yellow color fled color space. This table fills in entries within the color 

produced at the edge. The stripe of yellow color is not gamut of the target printer with the closest printable color, 

readily visible, however, when its luminance matches its 65 Colors out of the printer's gamut are represented by an 

green neighbor. However, as illustrated in FIG. 4, gamut averaging of the printable colors and of the outer range of 

mapping applied to the red/green edge will tend to decrease colors assigned for the saturates. A requested color is con- 
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verted to the reference color space, then looked up in the FIGS. 1, 2, 3 and 4 illustrate artifacts generated by out of 

table, and the resulting color specification is sent to the gamut processing; 

printer, FIG. 5 illustrates a printing system in which the present 

U.S. Pat. No. 5,450,216 to Kasson describes gamut- invention finds advantage; 

mapping color images in device-independent form to 5 FIG fi flustrates a functional block diagram illustrating 

device^dependent gamut in a perceptual color space. Digital one embodiment of the present invention, which may be 

images are mapped to any device-dependent gamut in a used m conjunction with the printing system of FIG. 5. 

manner that minimizes the human visual response both to F , G ? illustrates a Mft hue preserving gamut mapping 

the luminance and the chrominance changes necessary to memod usefijl ^ a . technique; 

force out-of-gamut pixels into the specified device- io ° . . 

dependent gamut. A "neighborhood gamut mapping" tech- FIG - 8 ^usfrates a nearest point gamut mapping useful as 

nique considers the subjective visual effects of nearby pixels a mapping technique; 

on the mapping of each pixel. At low spatial frequencies, FIG. 9 illustrates a gamut mapping method moving the 

image luminance is biased toward the luminance in the remapped points towards the provides gamut mapping 

device-dependent gamut at which the greatest chroma mag- 15 towards printer black and white; and 

nitude is available for a fixed hue angle. Spatial filtering FIG. 10 is a flow-chart illustrating the operation of a linear 

exploits the differing spatial frequency regions of insensitive adaptive spatial filter useful in conjunction with the inven- 

human visual response to both luminance and chrominance rion, 

changes. Referring now to the drawings where the showings are for 

The references cited herein are incorporated by reference 20 the purpose of illustrating an embodiment of the invention 

for their teachings. and not for limiting same, a basic system for carrying out the 

eiTx^AiDv /"it? roc i\T\/r\mnM present invention is shown in FIG. 5. In such a system, a 

SUMMARY OF THE INVENTION Source of original images 10, perhaps a color scanner, a 

In accordance with the present invention there is provided personal computer or workstation with appropriate docu- 

a method of gamut mapping that preserves variations in 25 ment or image creation software, a camera or data storage 

color within a local neighborhood of pixels in the image. device, produces an image, typically in a device independent 

In accordance with one aspect of the invention, there is manner provides suitable electronic images. The images are 

provided a method of remapping out-of-gamut colors in an defined in pixels, each associated with a color defined in 

original document, preparatory to use in a system having a terms of CIE color space L*,a*,b*, or some other 

different system gamut including: receiving a set of pixels luminance -chrominance space (L, C1,C2) or an analytic 

defining at least a portion of an image, each pixel mapped to transform thereof. 

an original pixel color; for any pixel in said set of pixels, An image from image da ta source 10 is directed to an 
remapping said pixel from said original color to a color ima ge processing unit (IPUT30, wEichTTor the purposes of 
which is within said system gamut; comparing colors in said 35 description includes all of the digital processing require- 
set of pixels before remapping with colors of said set of ments to convert an image from its original data represen- 
pixels after remapping and creating a comparison metric tation to a format suitable for printing at printer 30. Depend- 
therefrom; and applying said comparison metric to said ing on the image representation, image decomposition from 
colors after remapping to provide corrective feedback to said a PDL, halftoning from gray level images, color transforms, 
pixel set for final color selection within said system gamut. 4Q scaling and any other process required, is represented by the 

In accordance with another aspect of the invention, the IPU 20. IPU 20 can take various forms and operational 
comparison metric is filtered with a spatial filter prior to said details, and ranges from dedicated, and hardwired or soft- 
applying step. The spatial filter is selected to strengthen the ware driven operation in conjunction with high speed print- 
comparison metric in high frequency image regions to ing devices, to printer drivers running on personal comput- 
increase its impact on the gamut remapped colors, and 45 ers or workstations driving personal printers. Called out 
weaken the comparison metric in low frequency areas, to particularly is the gamut mapping function 24, which will be 
weaken its impact on the gamut remapped colors. further described below. It will no doubt be appreciated that 

It is an important aspect of the present invention that the me colors ma X be derived ™ r S b «> lor s P ace > md readn y 

mapping from an input color to an output color is not a converted to other device independent color spaces where 
unique mapping as done in standard gamut mapping tech- 5 o calculations can more readil y be performed. This may 

niques. Rather the mapping of an individual color depends readil y be part of the functionality of IPU 20. 

on its local neighborhood. This means that a specific color Printer 30 can be any color printer, printing a plurality of 

will be mapped to different output colors as a function of separations which, when superimposed, form a multicolor 

neighboring colors. image. For the purposes of the present invention and its 

The present invention seeks to preserve spatially local 55 description, the printer functionality could be replaced or 

luminance variations in the original image, by calculating arranged in parallel with a display or monitor, 

the difference between the original image luminance and the Gamut mapping function 24 serves to map pixels defined 

gamut mapped luminance resulting from using a pixel-bases by colors not printable by the printer, or resulting from initial 

gamut mapping process. A spatial filter is then applied to the image processing, into colors printable by the real printer. In 
difference signal, and the signal is added back to the 60 doin g *°> out-of-gamut, unprintable colors are mapped to 

gamut-mapped signal. The result is a reduction of the printable colors in accordance with some scheme that 

artifacts that arise from standard gamut mapping algorithms. attempts to optimize retention of the color relationships 

The invention may take form in various components and within the document. Colors that fall within the output 

arrangements of components, and in various steps and gamut may also be adjusted to retain the relationships to the 
arrangements of steps. The drawings are only for purposes 65 mapped colors. 

of illustrating a preferred embodiment and are not to be When referring to colorimetric or device independent 

construed as limiting the invention: spaces, the reference is to color space definitions that are 
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transforms of CIE XYZ space (1931). When we refer to 
device dependent space, we refer to a color space that is 
defined only in terms of operation of the device using it. 
While many color spaces have three dimensions, it is 
possible to have color spaces with less than three dimensions 
or more than three dimensions and likewise, it is possible for 
printers to use less than three colorants or more than four 
colorants. 

In the following discussion, the term "luminance" is used 
generically to encompass the strict definitions of luminance 
(i.e., the Y component in XYZ) and lightness (i.e. the L* 
component in L*a*b*). Chrominance components CI and 
C2 are likewise generic manifestations of opponent color 
signals, or the equivalent polar representations, chroma and 
hue. It should understood that exact input luminance values 
cannot always be preserved through the gamut mapping, 
since this can result in undesirable trade-offs in the other 
attributes, hue and chroma. The optimum trade-ofls in these 
attributes for various combinations of device gamuts, image 
types, and rendering intents is not the focus of the invention, 
and not discussed in detail here. 

Initially, it will be assumed that a gamut mapping function 
Gl is a reasonable pixelwise strategy for large color patches 
or smooth regions of an image. The invention preserves the 
characteristics of Gl at low spatial frequencies, while pre- 
serving high frequency luminance variations in the original 
image that Gl may have lost or altered. The second gamut 
mapping G2 is designed to approximately preserve these 
luminance variations. With reference now to FIG. 6, an 
image defined, for example in Y (luminance), CI, C2 
(chrominance space) is originally received and analyzed at 
gamut mapping function 100, which applies function Gl 
thereto. 

As noted, function Gl can be any gamut mapping 
function, perhaps represented by either an optimized non- 
linear function, or a simplistic mapping function, where out 
of gamut pixels are moved to the nearest "in gamut" plane. 
There are many such functions, and the precise function 
chosen is not relevant to the invention. The output from the 
function could be characterized as Y f , CI', C2\ From this 
output, the luminance signal Y' is compared to Y, perhaps in 
a simple differencing function, 102, producing AY In 
general, the other color channels could also be included in 
the differencing operation. A spatial filter 104 operating on 
an n x m block of values AY is used so that the feedback to 
the gamut mapping function is based on a local area, rather 
than derived on a pixel by pixel basis. The output of spatial 
filter 104, AY', is then used to vary Y* to obtain Y". In this 
example, AY is added to Y' at adder 106. The colors defined 
by Y", CI', C2' are then subjected to a second gamut 
mapping function 110, implementing gamut mapping func- 
tion G2, which may or may not be the same function as Gl. 
The resulting signals are directed on to printer correction 
function 112, for conversion to device dependent signals 
preparatory to printing. Gl and G2 may be provided as 
described in R. S. Gentile, E. Walowit, and J. R Allebach, A 
Comparison of Techniques for Color Gamut Mismatch 
Compensation, SPSE/SPIE Symposium on Electronic 
Imaging, Los Angeles, pp. 176-181; and E. G. Pariser,A/i 
Investigation of Color Gamut Reduction Techniques, IS&T 
2nd Symposium on Electronic Publishing, pp. 105-107, 
1991. 

Considering possible gamut mapping functions, several 
are possible and useful in the present application. The design 
of gamut mapping functions Gl and G2, and the spatial filter 
104 can depend on many factors, including global and local 
image characteristics, device characteristics, rendering 
intent and preference. 
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In one possible embodiment, gamut mapping function Gl 
was chosen to map out-of-gamut colors to the nearest 
surface point of the same hue, and G2 was chosen to map 
out-of-gamut colors to a surface point defined by the four 

5 steps: I) determine a hue of given out-of-gamut color; ii) 
calculate cusp point defined as the in-gamut color of maxi- 
mum chroma at the given hue; iii) determine focal point by 
projecting the cusp point onto neutral axis while preserving 
luminance; and iii) mapping the out-of-gamut color to 

10 surface in a direction towards the focal point. In this 
embodiment, for both Gl and G2, in-gamut colors were 
unaltered. The filter F was chosen to be a linear filter with 
support over an NxN image block, designed so that at low 
frequencies, AY'=0, while at high frequencies, AY'=AY This 

is is essentially a high pass filter. With these characteristics, the 
inventive gamut mapping process will approximately repro- 
duce the variations in Y at high spatial frequencies, while 
reducing to the pixelwise mapping Gl at low spatial fre- 
quencies. 

20 The examples illustrate where it is desirable to change the 
colors predicted by traditional gamut mapping in order to 
preserve information from the original image. In each case, 
the output color or colors selected by the proposed algorithm 
would have been different if the colors were large areas of 

25 solid colors. The desire in color reproduction of images is 
not to reproduce the closest colors to the original, but to 
reproduce the information and intent of the original. 
Other gamut mapping methods include: 

3Q With reference to FIG. 7, for a given input color, the hue 
angle thereof is obtained, and then the input color is mapped 
to the nearest point on the gamut surface within a range of 
hues about the input hue angle, where the range is greater 
than 0 and less than 90 degrees. A given point P, with hue 
angle H, is mapped to the closest point on the gamut surface 
within the region bounded by the lines W-Hl and W-H2. 

With reference to FIG. 8, a reduced gamut is defined by 
reducing the chroma of each surface point by a predeter- 
mined scale factor. For a given input color, a gamut-mapping 

40 vector is defined as the shortest distance to the reduced 
gamut. Map the input color to the original gamut surface in 
the direction given by this vector. 

With reference to FIG. 9, Point C is defined as the point 
of maximum chroma in the hue plane of the input color. 

45 Drawing two lines, WC from white to C, and KC from black 
to C. All points in the shaded region (e.g. P3) are mapped to 
point C. All points outside the shaded region whose lightness 
is greater than that of C (e.g. PI), are mapped to the surface 
in a direction towards black K. All points outside the shaded 

50 region whose lightness is less than that of C (e.g., P2), are 
mapped to the surface in a direction towards white W. 

Again referring to FIG. 6, the process applies pixelwise 
gamut mapping Gl to the input colors, and computes a 
difference image AY between the luminances of the input 

55 and gamut mapped signals. It then adds a filtered version 
AY of this difference back to the gamut mapped signal Y. 
This operation might move pixels back out of the gamut, and 
hence, a second gamut mapping operation G2 is needed to 
ensure that all colors are within gamut. Note that Gl and G2 

60 need not be the same algorithm. 

With reference now to FIG. 6 and spatial filter 104, one 
possible filter that finds use with respect to the invention is 
an adaptive linear filter. Using the above gamut mapping 
process, it can be seen that pictorial images require a large 

65 filter footprint in order to have a visual effect. Here "large" 
refers to roughly 15 to 30 pixels (1 -dimensional) for a 300 
dpi (dot per inch) printer. Smaller filters are not as effective. 
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Note that in general, it may be desirable to choose F as a 
nonlinear filter. The simple linear filter was chosen for ease 
of demonstration. It should also be noted that the filter can 
optionally have a gain larger than 1 to enhance detail. 

On graphical data, however, large filter sizes lead to a 
distinct and disturbing halo effect around edges. This effect 
is especially disturbing for text on a colored background. 
Consequently small filter windows or footprints should be 
used. When the wrong filter size is used, image quality 
improvement may not be obtained. 

The filter can be adapted to improve its use by changing 
the filter footprint as a function of the local image data. In 
this way, the filter function adapts to local image content. It 
should be noted that one could also change the filter values 
or coefficients (rather than the footprint). However, variation 
of the filter footprint is a preferred option, because a fast 
implementation exists for a simple, constant-weight blur 
filter. 

One way of looking at the filter selection issue can be 
described in the framework of pictorial vs. graphic data. A 
closer examination, however, reveals that the issue is con- 
nected to the local image data variation in terms of local 
activity and "noise". Clean edges in the input image data are 
best mapped using a smaller filter footprint. Noisy and weak 
edges are best mapped using a large filter footprint. A 
significant improvement is obtained by using different filter 
footprints based on local image characteristic. Different 
local image activity measures can be used, but they all have 
in common that: 

a) Very low activity [flat areas] correspond to a small filter 
footprint (=1; 5x5 footprint, for example). 

b) Medium activity [mild image variations] correspond to 
a large filter footprint (15x15 footprint, for example). 

c) High activity [strong edges] correspond to a small filter 
footprint (>1; 5x5 footprint, for example). 

The most obvious metric is data norm of order p: 

Here, e,- is the luminance error at the center pixel i and e ; - is 
the luminance error at pixel j within the predefined neigh- 
borhood of pixel i. The most common forms of this activity 
measure are: 
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For speed of implementation, we are using a modified 
activity measure: 



55 



The advantage of the above activity measure lies in the fact 
that the quantity used as activity is already computed as part 
of the standard filtering described. 

The described filtering function may be implemented by 
generating an activity metric at activity metric function 400, 
which applies its output to filter selector 302, to select an 
appropriate filter for spatial filter 104. 

FIG. 7 shows a flowchart of the adaptive filtering method. 
At step 300, a small area activity metric as is computed, 
conveniently using the absolute value of the small area filter 
calculation. At step 302, if the small area activity a, falls 



60 
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outside a medium range, the small area filter is used (step 
304). If & s falls within the range, a large area activity z L is 
computed at step 306, again conveniently using the absolute 
value of the large area filter calculation. At step 308, values 
of the activity metrics are compared. If the ratio of the two 
activity metrics a/a^ is larger than a specified threshold T 3 , 
a small filter is used at step 310; otherwise a large filter is 
used at step 312. It should be noted that re -normalization 
constants in the activity metrics based on the filter size have 
been omitted. 

The required filters vary with the object or image type of 
the input, but simple activity measures as described above 
are sufficient to select the appropriate filters. It should be 
clear that the selection between two discrete filters was 
chosen for simplicity and speed and that in the general case 
the activity measure can be used to derive a filter or to 
smoothly blend between a number of pre-selected filters. 

It Will no doubt be appreciated that the present invention 
can be accomplished through application software accom- 
plishing the functions described, to operate a digital com- 
puter or microprocessor, though a hardware circuit, which 
will probably provide optimum speed, or though some 
combination of software and hardware. 

The invention has been described with reference to the 
preferred embodiment. Obviously, modifications and alter- 
ations will occur to others upon reading and understanding 
the preceding detailed description. It is intended that the 
invention be construed as including all such modifications 
and alterations insofar as they come within the scope of the 
appended claims or the equivalents thereof. 

What is claimed is: 

1. A method o remapping colors in an original document, 
preparatory to use in a system having a different system 
gamut, comprising the steps of: 

receiving a set of pixels deferring at least a portion of an 
Image, each pixel associated with an original pixel 
color; 

for any pixel in said set of pixels, remapping said pixel 
from said original pixel color to a new pixel color 
within said system gamut; 

comparing original pixel color in said set of pixels with 
new pixel colors in said set of pixels and creating a 
Comparison metric therefrom, wherein said compari- 
son metric is a difference between an original pixel 
color luminance component and a new pixel color 
luminance; 

processing said comparison metric through a spatial fil- 
tering operation; and 

combining an output of said spatial filtering operation 
with said new pixel colors. 

2. A method as described in claim 1, wherein said remap- 
ping is accomplished by remapping each pixel from said 
original pixel color to said new pixel color, said new pixel 
color emphasizing preservation of chroma or saturation over 
preservation of lightness of the pixel. 

3. A method as Described in claim 1 where said combin- 
ing adds an output of said spatial filtering to said new pixel 
color luminance component. 

4. A method as described in claim 1 where said spatial 
filtering uses a linear filter. 

5. A method as described in claim 1 wherein said spatial 
filtering strengthens said comparison metric in high spatial 
frequency regions to increase the effect thereof on new pixel 
colors, and weakens said comparison metric in low spatial 
frequency regions to weaken the effect thereof on new pixel 
colors. 
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6. A method as described in claim 1, further comprising an 
additional remapping of said set of pixels associated with 
new pixel colors to adjusted new pixel colors, where said 
combining has moved said new pixel colors out of said 
system gamut. 

7. A method as described in claim 6, wherein said addi- 
tional remapping comprises: 

determining a hue of each said new pixel color; 
calculating a cusp point, defined as a point on a printer 

gamut surface corresponding to a maximum achievable 

chroma for said hue; 
determining a focal point, defined as a projection of said 

cusp point into a neutral axis of a color space, in a 

manner that preserves luminance; and 
mapping pixel colors that are out of said system gamut to 

said adjusted new color on a surface of said system 

gamut in a direction in said color space towards said 

focal point. 

8. A method of remapping colors in an original document, 
preparatory to use in la system having a different system 
gamut than originally applied, including: 

receiving a set of pixels defining at least a portion of an 
image, each pixel therein mapped to a first set of 
original colors. 

for each pixel in said set, remapping said pixel from said 
set of original colors to a new set of colors within said 
system gamut; 

determining a comparison metric, comparing said set of 
pixels mapped to said first set of original colors and 
said set of pixels mapped to said new set of colors, 
wherein said comparison metric is a difference between 
luminance values of said set of pixels mapped to said 
first set of original colors and said set of pixels mapped 
to said new set of colors; and 

applying said comparison metric to said colors after 
remapping to provide corrective spatial feedback to 
said pixel set to adjust color selection within said 
system gamut. 

9. The method as defined in claim 8, and including 
filtering said comparison metric prior to corrective feedback. 

10. The method as defined in claim 8, wherein said 
filtering is provided by a linear filter. 

U. The method as defined in claim 8, wherein said linear 
filter strengthens said comparison metric in high spatial 
frequency regions, to increase the effect thereof on said new 
set of colors, and weakens said comparison metric in low 
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spatial frequency regions, to weaken the effect thereof on 
new set of colors. 

12. The method as defined in claim 8, wherein said 
filtering is provided by a non-linear filter. 

13. The method as defined in claim 8, and including an 
additional remapping of said pixels to said new set of colors, 
to correct any mappings where said corrective feedback has 
moved said new color out of said system gamut. 

14. A gamut mapping system, including: 

an image: input, receiving pixels directed thereto mapped 
to a first set of colors, which may not be reproducible 
at a selected output device; and 

a first gamut mapping device, remapping said pixels 
mapped to said first set of colors to a second set of 
colors reproducible at said selected output device; 

a corrective spatial feedback system, having: 

a first input receiving said pixels mapped to said first set 
of colors, 

a second input receiving said pixels mapped to said 
second set of colors, and producing a comparison 
signal representing a difference therebetween, said 
difference being between an original pixel color 
luminance component and a new pixel color lumi- 
nance; and 

a signal ladder, adding said comparison signal to a lumi- 
nance value of said pixels mapped to said second set of 
colors, to adjust said luminance values, while main- 
taining other color defining values. 

15. A gamut mapping system as defined in claim 14, and 
including: 

a spatial filter, operating on said comparison signal before 
it added to said luminance value of said pixels mapped 
to said second set of colors, said spatial filter strength- 
ens said comparison metric in high spatial frequency 
regions, to increase the effect thereof on said new set of 
colors, and weakens said comparison metric in low 
spatial frequency regions, to weaken the effect thereof 
on new set of colors. 

16. A gamut mapping system as defined in claim 14, and 
including 

a second gamut mapping device, remapping said pixels 
mapped to said second set of colors, after said com- 
parison signal is added thereto, to assure said pixels 
mapped to said second set of colors are reproducible at 
said selected output device. 
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